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ABSTRACT 

IGR J04571+4527 and Swift J0525.6+2416 are two hard X-ray sources detected in the 
5wt/f/BAT and INTEGRAL/IBIS surveys. They were proposed to be magnetic cataclysmic 
variables of the Intermediate Polar (IP) type, based on optical spectroscopy. IGR 104571+4527 
also showed a 1218 s optical periodicity, suggestive of the rotational period of a white dwarf, 
further pointing towards an IP classification. We here present detailed X-ray (0.3-10 keV) 
timing and spectral analysis performed with XMM-Newton, complemented with hard X-ray 
coverage (15-70 keV) from 5w//f/BAT. These are the first high S/N observations in the soft 
X-ray domain for both sources, allowing us to identify the white dwarf X-ray spin period 
of Swift J0525.6+2416 (226.28 s), and IGR 104571+4527 (1222.6 s). A model consisting of 
multi-temperature optically thin emission with complex absorption adequately fits the broad¬ 
band spectrum of both sources. We estimate a white dwarf mass of about 1.1 and 1.0 Mq for 
IGR 104571+4527 and Swift 10525.6+2416, respectively. The above characteristics allow us 
to unambiguously classify both sources as IPs, confirming the high incidence of this subclass 
among hard X-ray emitting Cataclysmic Variables. 

Key words: novae, cataclysmic variables - white dwarfs - X-rays: individ¬ 

ual: Swift 10525.6+2416 (also known as IRXS 1052523.2+241331 and as WISE 
1052522.84+24133333.6), IGR 104571+4527 (also known as IRXS 1045707.4+452751 
and as Swift 10457.1+4528). 


1 INTRODUCTION 

Binary systems containing an accreting white dwarf (WD) from a 
late-type main sequence or sub-giant star are known as Cataclysmic 
Variable stars (CVs). Thanks to the recent hard X-ray surveys car - 
ried out by Swift Burst Alert Telesc ope (B AT; Ib arthelmTbOOOl) 
and INTEGRAL IBIS dUbertini et al.ll2003h instruments above 20 
keV, the number of det ected magnetic sys t ems (where B > 10^ 
rrlow et al.l2006l:[Bemardini et alj2oT^ 

therein). Magnetic CVs are about 20 per 
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G) is rapidl y increasing dB; 
|2013L|2 oT^ and references 


cent of all CVs dPerrario. de Martino & Gansickell20f^ and they 
are further divided in two main groups, depending on the degree 
of synchronization (Prot^ui/Porb=n]- The synchronous systems, 
called polars, are polarised in the optical and near-IR bands re¬ 
vealing WDs with strong magnetic fields (B ~ 7 — 230 x 10® 
G). These are able to lock the rotation of the WD with the orbital 
period. Those systems which instead possess an asynchronously 
rotating WD are called intermediate polars (IPs) and generally do 
not show optical/nearIR polarisation, suggesting they are weakly 
magnetised WDs (B < 10® G). In magnetic CVs, the accre¬ 
tion flow is magnetically channelled on the WD polar regions. 
The high field Polars accrete matter directly through a stream 
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from the Roche lobe overflowing donor star, while the lower 
field IPs may accrete through a stream or a truncated disc or 
ring, dependin g on their degree of asynchronism a n d their mag¬ 
netic moment dNorton, Wynn & SomerscalesI l2004l : iNorton et alJ 
llOOSh . The the magnetically confined accretion flow approaches 
supersonic velocities and a strong stand-off shock 10 — 80 
keV) is formed. In the post-shock region matter slows and cools 
via cyclotron radiation emerging in the optical/ne ar-IR and/or 
via thermal bremsstrahlung emerging in the X-ravs dAiz t3 ll973l: 


_ tug I __ 

IWu. Chanmugam & Shavi^ll994l : ICroDper et alJll999h . The dom¬ 
inant cooling channel primarily depends on the WD magnetic 
field intensity dWoelk & Beuermannlll99^ : Ipischer & BeuermannI 
l200lh . The post-shock region in the polars mainly cool via cy¬ 
clotron emission, while IPs are bremmstrhalung dominated sys¬ 
tems. The post-shock emission is intercepted and thermalised by 
the WD surface, and can appear to a far observer as an optically 
thick, highly absorbed, spectral component emerging in the soft 
X-ray (~ 20 — 60 eV) or in the ultraviolet or extreme ultraviolet 
(UV/EUV). 

Magnetic CVs are debated as important contributors of the 
X-ray luminosity function at low luminosities (10®° — 10®® 
erg/s) as inferred from X-ray surveys of the galactic ridge and 
bulge I Sazonov et al. 20061: ^vnivtsev et alj l2008l . I 2 OO 9 II 20 III : 
iHong et alj I 2 OI 2 I : iHond 1201 A . The fact that they are about 
6% of the galactic hard X-ray sources in th e Swift/BAT and 


INTEGRAL/IBIS sur veys dBarlow et al. 
ICusumano et alj|20ld : [^umgartner et al 


20061: iBirdetaP I 2 OI 0 I: 
supports this sce- 
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narios. IPs represent about 80% of these hard X-ray CVs, and 
their number is increasi ng thanks to followmps in the X-rays and 
optical bands (s ee e. g. iBarlow et alj I 2 OO 6 I : lAnzolin et aljl200^ : 
[Bernardini et alj|2012ll2013 1. The characterisation of new sources, 
especially the faint ones, may allow inferring new aspects of 
the luminosity function, such as a possible large, but still to-be- 
discove red, population of low-lu minosity and short orbital period 
IPs (see lPretorius & Mukaill201^ . It will also allow understanding 
the role played by main source parameters, such as the WD mass, 
magnetic field intensity, and the mass accretion rate in shaping the 
evolution of these systems and their X-ray spectral properties. X- 
rays are the best diagnostic bands to identify the true WD spin pe¬ 
riod, thus allowing a secure classification. Moreover, new proper¬ 
ties, such as the presence of a soft blackbody component and of a 
ionised (warm) absorber have been recently identified in an increas¬ 
ing number of IPs thank s to the unique capabilities of the XMM- 
Newton satellite (see e .g. lMukai et alj200ll : lde Martino et alj2008l : 
[Bernardini et alj2012t) . 

We here present the first high S/N soft X-ray observations, car¬ 
ried out with XMM-Newton, of Swift J0525.6-1-2416 (also known as 
IRXS J052523.2-H241331 and as WISE J052522.84-H24133333.6) 
and IGRJ04571-1-4527 (also known as IRXS J045707.4-H452751 
and as Swift J0457.1-1-4528). Based on the characteristics of its 
optical spectra. SwiftJ0525.6-1-2416 was prop osed to be a mag¬ 
netic CV or a high accretion rat e novalike jTorres et alj 1 20071 : 
iMasetti et alj2013) . lRams^ ( l2009h reported the possible detection 
of an eclipse from this source by using optical spectroscopic data. 
IGR J04571-t4527, inste ad, was first classified as a non-magnetic 
CV iMasetti et alj|2010h . Later on, the unambiguous detection of 
a coherent signal at 1218.7 ± 0.5 s in optical photometric data, 
interpreted as the WD spin period, together with the detection of 
a spectroscopic period (4.8 h or 6.2 h, one the alias of the other: 
iThorstensen & Halpeml[2oI^ . provided strong arguments on the 
IP nature of IGR J04571-H4527. 

We complement the X-ray study of both sources with 


simultane ous optical photo metry acquired with the Optical 
Monitor iMason et ^ [ 20 ^ on board XMM-Newton and for 
Swift J0525.6-1-2416 we also i nclude a re-ana lysis of the optical 
spectroscopic data presented in iRamsa 3 i l2009l) and Wide-field In¬ 
frared Survey Explorer (WISE) archival photometry. 


2 OBSERVATIONS AND DATA ANALYSIS 
2.1 XMM-Newton observations 

SwiftJ0525.6-l-2416 and IGRJ04571-1-4527 were observed on 21 
February 2014 and 13 Marc h 2014, respectively , with the EPIC PN, 
MOST and MOS2 cameras jStriider et al.l20oi] : lTumer et al.l200lll 
as main instruments. The log of the observations is reported in Ta- 
ble[T] Data were processed using the SAS version 13.0.0 and the 
latest calibration files (CCF) available on October 2014. 


2.1.1 The EPIC data 

We extracted the source photons from a circular region with a ra¬ 
dius of 40 arcsec, centered at the source position. Background pho¬ 
tons were taken from a nearby empty region of the sky, with a ra¬ 
dius of 70 (PN) and 140 (MOSs) arcsec. For Swift J0525.6-I-2416 
we used the whole data set for both the spectral and the timing 
analysis (the high background contamination is only a small frac¬ 
tion of the source flux). In the case of IGRJ04571-1-4527, the high 
background epoch is mainly concentrated at the end of the point¬ 
ing, consequently we removed this part of the observation. For the 
timing analysis, instead, we used the whole data set. 

We produced background-subtracted light curves in the ranges 
0.3-12 keV, 0.3-2 keV, 2-3 keV, 3 - 5 keV and 5 - 12 keV (for 
Swift J0525.6-1-2416 we also used the 0.3 — 1 and 1 — 3 keV range). 
The event arrival times were corrected to Solar system barycentre 
by using the task BARYCEN. To study spectral variability with the 
source rotation we extracted the EPIC spectra at pulse maximum 
and minimum. 

The EPIC spectra were rebinned before fitting with the tool 
SPECGROUP in order to have a minimum of 50 (PN) and 25 (MOS) 
counts each bin. The spectral fit were made with XSPEC version 
12.8.2. The sources are too faint for a use of the RGS spectra. 


2. 1.2 The Optical Monitor photometry 

We run the task OMECHAIN to produce the B-band light 
curves. However, due to offset coordinates, the pointing to 
IGRJ04571-I-4527 had the fast window centred on a different star. 
Thus, no fast timing on this source is available, but magnitudes 
were obtained from the 2-D images that are acquired along each of 
the ten ~3200 s long segments. In these images IGRJ04571-1-4527 
is stable at an average magnitude B= 18.35 ± 0.03 corresponding 
to a flux of 3 X 10“®® erg c m~® s~® A~®. It is 1.6 times fainter than 
when observed in 2009 by IMasetti et alj 1 2010h and at about the 
same l evel than when observed in 2010 by Thorstensen & Halpernl 
( I 2 OI 3 I) . Swift J0525.6-1-2416 has an average magnitude B= 17.57± 
0.01 corresponding to a B-band flux of 6.2 x 10“®® erg cm~^ s~® 
A“®, which is comparable with what found bv iMasetti et all 1 I 2 OI 2 I) 
and reported in the USNO A2.0 catalogue. 
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Table 1. Summary of main obsei'vations parameters for all instruments. Uncertainties are at Icr confidence level. 


Source 

Telescope 

OBSID 

Instrument 

Date 

yyyy-mm-dd 

UTstart 

hh:mm 

'-ra 

^ expo 

(ks) 

Net Source Count Rate 

c/s 

IGR104571+4527 

XMM-Newton 

0721790201 

EPIC-PN^ 

2014-03-13 

11:42 

36.1 

1.424 ±0.008 




EPIC-MOSL 

2014-03-13 

11:19 

37.8 

0.491 ± 0.004 




EPIC-MOS2= 

2014-03-13 

11:19 

37.8 

0.494 ± 0.004 




OM-B'^ 

2014-03-13 

11:25 

32.0 

2.33 ± 0.07 


Swift 

e 

BAT 



8.5 X 10® 

3.9 ±0.3 X 10-® 

SwiftJ0525.6+2416 

XMM-Newton 

0721790301 

EPIC-PN^ 

2014-02-21 

07:15 

30.0 

1.494 ±0.008 




EPIC-MOSL 

2014-02-21 

06:52 

31.7 

0.574 ± 0.04 




EPIC-MOS2= 

2014-02-21 

06:52 

31.7 

0.577 ±0.04 




OM-B 

2014-02-21 

06:58 

26.6 

4.78±0.03 


Swift 

e 

BAT 



7.3 X 10® 

3.7 ±0.4 X 10-5 


“ Net exposure times. 

^ Full frame mode (thin filter applied). 

Large window mode (thin filter applied). 

Fast window mode. The central wavelength of the filter is 4500 A. 

® All available pointings collected during 2004 December to 2010 September are summed together. 


2.1.3 The optical spectroscopy of Swift J0525.6+2416 

Spectroscopic observations of SWIFT J0525+2413 were obtained 
using the 4.2m William Herschel Telescope (WHT) equipped with 
the Intermediate dispersion Spectrograph and Imaging System 
(ISIS) on La Palma. Four spectra, each with an exposure time of 
240 sec, were taken consecutively on the night of 2008 Oct 6. The 
cadence was ~9 min with a start-to-end duration of ~32 min. Con¬ 
ditions on the night were clear, hut with variable seeing. We used 
the R300B (~3500-5300 A) and R158R (~5700-10000 A) grat¬ 
ings together with a 0.8 arcsec slit. A fit to the arc lines gave a =1.9 
and a =2.2 A for the blue and red spectra, respectively. The data 
were reduced using the StarlinkQ Figaro package. The wave¬ 
length calibration was obtained using CuAr-l-CuNe arc lamp expo¬ 
sures taken before the series of spectra on the target. Observations 
of the spectroscopic standard He3 were used to remove the instru¬ 
mental response. Due to the variable seeing during the night some 
flux could have been lost and therefore we did not calibrate the 
absolute flux scale. 


2.2 The Swift observations 

Both sources are catalogued in the 70 month Swift/BAT 
source catalogue (Baumgartner et al. 2013). We down¬ 
loaded the archival eight-channel spectra of the two sources 
and the response matrix files that are directly available at 
http://swift.gsfc.nasa.gov/results/bs70mon/ Due to the low S/N 
of the BAT spectra above 70k eV we limited the analysis of the 
broad-band spectra up to this energy. 


^ The Starlink Software Group homepage can be found at 
http://starlink.jach.hawaii.edu/starlink 


3 RESULTS 

3.1 IGRJ04571-I-4527 

3.1.1 The X-ray and optical timing analysis 

The 0.3-12 keV PN background subtracted light curve shows a fac¬ 
tor of ~ 2.5 variability on timescales of tens of minutes through all 
the pointing. There is no evidence for long term variability asso¬ 
ciated with the optical 5-6 h orbital period. The power spectrum 
of the combined PN plus MOSs 0.3-12 keV light curve shows a 
peak at ~ 8 X 10“^ Hz (see Fig.[T]l. The exact period of the main 
peak was determined by means of a phase-fitting technique (see 
iDall’Osso et al.ll2003l . for more details on this technique), and it is 

= 1222.6 ± 2.7 s. All uncertainties are hereafter at la con¬ 
fidence level. This is consistent within 1.5 sigma with the optical 
period reported by iThorstensen & Halpeml l l2013l) . We do not de¬ 
tect significant power at higher harmonics. To further test that the 
period is the true spin period of the WD, we folded the light curve 
at twice the 1222.6 value without finding changes in the pulse in 
two consecutive cycles. We conclude that 1222.6 is the period of 
the fundamental harmonic and thus it is the WD spin period. 

We also searched for variations of the phase of the spin pe¬ 
riod with time and found a sinusoidal trend (see Fig.[TJ. We infer 
a period of T.QiJ;® h, consistent within 2a with the longest (6. 2 
h) spectroscopic period found by IThorstensen & Halpernl ( l2013h . 
We stress that this period is consistent with the total duration of the 
pointing and that one phase-bin does not lie on the sinusoidal trend 
(we did not include it in the fit). 

We also studied the dependence of the pulses with the 
energy interval by folding the background subtracted light 
curves at the more accurate period found by period found by 
IThorstensen & Halpeml ^20131) , 1218.7 ± 0.5 s. A sinusoidal fit to 
the 0.3-12 keV folded light curve gives a pulsed fraction (PF) of 
10.0 ± 0.5 per cent. Here the PF is defined as: PF = {Amax — 
Amin)/{Amax + Amin), where Amax and Amin are the maxi¬ 
mum and minimum value of the sinusoid, respectively. In Fig. [T] we 
show the 0.3-12 keV pulse profile. We found that the PF slightly 
decreases with the energy interval: PFo.3-2 = 12.0 ±0.8 per cent, 
PF2-3 = 11.0 ± 1.2 per cent, PF3-3 = 9.8 ± 1.2 per cent, and 
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Figure 1. Up'. PN plus MOSl and MOS2 0.3-12 keV power spectrum of 
IGRJ04571+4527. The peak at ~ 8 x 10“^ Hz corresponding to the WD 
spin period is highlighted with an aiTow. Center. Evolution of the phase of 
the signal corresponding to the WD spin period versus time. The measured 
sinusoidal period is also reported in the upper right comer. The red points 
do not follow the sinusoidal trend and it is not included in the fit. Down'. 
X-ray pulse profile in the energy interval 0.3-12 keV. The folding period is 
1218.7 s. Two cycles are shown for plotting purposes. 


PF 5-12 = 7.3 ± 1.3 per cent. We also inspected the phase depen¬ 
dence of hardness ratios hut did not find any significant variations. 


3.1.2 The X-ray spectral analysis 

A single thermal model does not fit adequately the hroad-band 
(0.3-70 keV) average spectrum. Instead, we found a satisfying 
spectral fit by using a model composed by the sum of two MEKAL 
(multi-temperature optically thin plasma) and a Gaussian at 6.4 
keV, accounting for the fluorescent iron line. The spectral model 
also requires a complex absorption, made by a total PHABS plus a 
partial covering PCFABS absorptions, as frequently found in mag- 


Table 2. IGR104571+4527 spectral parameters for the best fitting model. 
We report the absorbed/(unabsorbed) 0.3 —lOkeV flux (Fq. 3-10), together 
with the 10 — 100 keV flux (Fio—loo)- Uncertainties are at the Icr confi¬ 
dence level. 


IGR J04571+4527 average spectrum 

Nfl-p 


10^^ cm~2 

0.302 ± 007 

Nffpc 


10^2 cm-2 

5.4 ± 1.0 

cvf 


% 

32 ± 3 

kThot 


keV 

63±f9 

nornihot 


10-3 

5.8 ±0.4 

kTcoid 


keV 

6.8±\'l 

nornicoid 


10-3 

1.3 ±0.4 

Az 



0.87 ±0.20 

EW“ 


keV 

0.13 ± 0.01 

Fo.3-10 

10“ 

-1 1 —9 —1 

erg cm s 

0.973 ± 0.010 (~ 1.3) 

FlO-lOO 

10“ 

-1 1 —9 —1 

erg cm s 

1.9 ±0.2 

xl (dof) 



1.10(362) 


“ Gaussian energy fixed at 6.4 keV. 


netic CVs (see e.g. lBemardini et ^l2012h . The use of a partial ab¬ 
sorber is dictated by the slight increase of PF at lower energies, 
suggestive of localised absorbing material. The results of the spec¬ 
tral fit are reported in Tab. and shown in Fig. The hydrogen 
column density of the total absorber is lower than the total galactic 
colum n density in the direction of the source (5.4 — 6.5 x 10^^ 
cm~^. loickev & Lockmanl[l990l : iKalberla et alj|2005h . indicating 
an interstellar origin. The Ah of the partial (32 percent) cover¬ 
ing absorber is about an order o magnitude larger, suggesting that 
it has a local origin. The spectral fit indicates that the post-shock 
flow is dominated by rather high temperat ures, without requiring 
a low-temperature plasma component (see [Bernardini et al.ll2012L 
and Sect. 3.2 of this work). We notice that a soft BB component is 
not required by the fits. We also fitted the broad-band spectrum wit h 
the model developed by Suleimanov, Revnivtsev & Ritterl (l2005h . 
which takes into account both temperature and gravity gradients 
within the post-shock region. This more physical model, imple¬ 
mented in Xspec, allows a more reliable estimate of the shock tem¬ 
perature and thus the WD mass. To avoid the effects of the com¬ 
plex absorption we fitted the spectrum above 3 keV only and added 
two Gaussians lines that take into account the fluorescent iron line 
and the thermal complex. We derive Mwd = 1.12 ± 0.06 M© 
(xS = 1.13, 215 d.o.f.). This value is slightly larger but consistent, 
within uncertainty, with that derived using the hot MEKAL temper¬ 
ature (1.03 M©). 

The spectra at pulse minimum (<j) = 0.12 — 0.37) and maxi¬ 
mum (<j) = 0.60 — 0.90) were fitted using using the spectral model 
reported in Tab|2 fixing the interstellar absorption and the abun¬ 
dance at their best fitting average value. We also fixed the hot 
MEKAL temperature to the average value, as it resulted to be un¬ 
constrained. This is not surprising given the weak PF at high en¬ 
ergies. We get xi ~ 1-01> 274 d.o.f. and xl = 1-03, 246 d.o.f. 
for the pulse maximum and minimum respectively. While all other 
parameters are constant within statistical uncertainties and consis¬ 
tent with the values found for the average spectrum, the only two 
parameters that change at maximum and minimum of the pulse are 
the Nh of the partial covering absorber and the normalisation of 
the cold MEKAL. The former is 10 ± 3 x 10^^ cm“^ at minimum, 
whilst it decreases to 4.7 ± 0.8 x 10^^ cm“^ at maximum. The lat- 
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Figure 2. Left panel: IGRJ04571+4527 broad-bad count spectrum. PN spectrum is in black, MOS 1 in red, MOS 2 in green and BAT in blue. Residual are 
shown in the lower panel. Right panel: The same as the left panel, but for Swift J0525.6+2416. 


ter instead is 2.0 ± 0.1 x 10 ® at pulse maximum, while decreases 
to 0.7 ± 0.3 X 10“® at minimum. 


3.2 Swift J0525.6+2416 

3.2.1 The X-ray and optical timing analysis 

We first visually inspected both the optical B band and the 0.3- 
12 keV PN background subtracted light curves using a small bin- 
ning sample (1 0 s) to check for sign of the eclipse reported in 
iRamsavl i2009li and we did not find any. While the X-ray light 
curve is showing short term variability, the B band one is almost 
flat with an average count rate of 4.8 c/s. The power spectrum of 
the combined PN plus MOSs 0.3-12 keV light curve shows a nar¬ 
row peak at ~ 4.5 x 10“® Hz (see Fig. central panel) with 
indication of a weaker harmonic at ~ 9 x 10“® Hz. The period 
of the main peak, determined with the phase-fitting technique, is 

_ 226.28 ± 0.07 s. This is the first detection of a coherent 
signal from this source. There is no sign of orbital variability in the 
evolution of the phase of the main peak with time. Also in this case 
we checked that the 226.3 s period is the fundamental. We also do 
not find any long term trend (> 300 s) in the X-ray light curve. 
Thus, we conclude that the spin period of the source is 226.3 s, and 
there is no evidence of eclipses. 

The PF of the 0.3-12 keV light curve for the fundamental and 
first harmonic is 6.3 ± 0.5 per cent and 3.1 ± 0.5 per cent, re¬ 
spectively. In Fig. Owe show the pulse profile for the 0.3-12 keV 
interval. The PF slightly decreases as a function of the energy: 
PFi -2 = 10.3 ± 1.1 per cent, PF 2-3 = 4.6 ± 1.6 per cent, 
PF 3-3 = 4.7 ± 1.4 per cent, and PF 3-12 = 5.4 ± 1.4 per cent. 
We did not include the 0.3-1 keV interval because of the low S/N. 
We also did not find any variation of the hardness ratios with phase. 


3.2.2 The X-ray spectral analysis 

The broad-band 0.3-70 keV spectrum requires more than one ther¬ 
mal component. We first used a model composed by the sum of two 
MEKAL and a Gaussian fixed at 6.4 keV absorbed by a total PHABS 
plus a partial covering PCFABS absorber. We found residuals (in 
excess with respect to the model) at low energy around 0.55-0.65 
keV, where OVlll and OVll lines are expected, respectively. We ob¬ 
tained a better fit by substituting the two MEKAL components with 
two VMEKAL, where the abundance of the different elements can be 
individually controlled. We left the oxygen abundance free to vary, 
linking the abundances of all other elements to a single, global, 
value and and found Ao = 0.14 ± 0.05 compared to the Solar 
value. The abundances of elements other than oxygen, were found 
to be Az = 0.52 ± 0.2(0. Summarizing, the MEKAL model pre¬ 
dicts an overabundance of oxygen, this problem is resolved using 
the VMEKAL model that suggests a subsolar oxygen abundance. We 
also notice that a soft BB component is not required by the fits. The 
yd/ is not satisfactory (1.29 for 395 d.o.f.), althoug h this can happen 
when dealing with PN and/or MOS data (see e.g. [Bernardini et alj 
l2012h for this kind of sources. We have tried several alternative 
models but none resulted in an improved yd- The results of the 
spectral fit are reported in Tab. and shown in Fig. We notice 
that the hydrogen column density of the total absorber is lower 
than the total galactic co lumn density in the direction of the source 
(3.2 — 3.8 X 10^^ cm^ [Pickev & Lockm^ll99d : iKalberla et alj 
l2005h . The tem perature of the cool ME KAL is similar to those found 
in other IPs (see lBernardini et alj2Q12h . The WD mass derived with 


^ The abundances we derive depend critically on the abundance ta- 
ble we adopt. However, subso lar ab undances are confirmed us ing both 
Iwilms. Allen & McCravl i2000l) and lAnders & Greyest jl989h relative 
abundances, being the former a factor of about two larger. Since the abun¬ 
dances of O and Fe, the two elements that our data are most sensitive to, 
differ greatly between the two tables, this is not surprising. 
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Frequency (Hz) 



Phase 


Figure 3. Up'. PN plus MOSl and MOS2 0.3—12 keV power spectrum of 
Swift J0525.6+2416. The nan'ow peak at ~ 4.4 X 10“® Hz, con'esponding 
to the WD spin period, is highlighted with an arrow. The first harmonic is 
also evident. Down'. 0.3-12 keV pulse profile. The folding period is 226.28 
s. Two cycles are shown for plotting purposes. 


the model presented in ISuleimanov. Revnivtsev & Ritteij (l2005l) is 
Mwd = 1-01 ± 0.06 Mq (xl = 1.18, 231 d.o.f.). On the other 
hand, we get Mwd = O.SSioiof M© when using the hot VMEKAL 
temperature. The lower values also reflect the use of only one- 
temperature (instead of a gradient of temperature) for the hottest 
post-shock regions. 

The spectrum at pulse minimum ((j> = 0.45 — 0.69) and maxi¬ 
mum (</> = 0.85—1.05) were fitted with the model presented in Tab. 
[3 fixing the interstellar absorption and the abundance at their best 
fitting average value. We also fixed the hot VMEKAL temperature 
to the average value, as it resulted to be unconstrained. We obtain 
~ 1.15, 282 d.o.f. and xS = 1.12, 289 d.o.f. for the pulse max¬ 
imum and minimum, respectively. We found that the the density 
of the localized absorber is slightly lower at pulse minimum than at 
maximum, 14.6±2.0 x 10^^ cm“^ with respect to 25.8±4.0 x 10^^ 
cm~^, and that the normalisation of the hot VMEKAL is slightly 
higher at pulse maximum than minimum, 0.0164 ± 0.0005 with 
respect to 0.0092 ± 0.0003. The cold VMEKAL parameters and the 
covering fraction are found consistent, within uncertainties, with 
their average values. 


3.2.3 Optical spectroscopy 

iRamsavl ) l2009l) reported that three of the four WHT spectra show 
B aimer series and Hell (4686 A) in emission. The second spectrum 


Table 3. Swift J0525.6+2416 spectral parameters for the best fitting model. 
We report the absorbed/(unabsorbed) 0.3 — 10 keV flux (Fq.s— lo) together 
with the 10 — 100 keV flux (Fio—too)- Uncertainties are at the la confi¬ 
dence level. 


Swift J0525.6+2416 average Spectrum 


N+fp 


10 ^^ cm“2 

0.31 ±0.02 

Nffpc 


10^2 cm-2 

18.1 ± 1.5 

cvf 


% 

45 ±2 

kThot 


keV 

40±i 

nomihot 


10-3 

10.0 ±0.7 

kTcold 


keV 

0.19 ± 0.01 

nonueoid 


10-3 

6.0±3;3 

Az“ 



0.52 ± 0.20 

Ao “ 



0.14 ±0.05 

EW*’ 


keV 

0.12 ± 0.03 

Fo.3-10 

10 - 

-1 1 —9 —1 

erg cm s 

1.13 ±0.02 (~ : 

FlO-lOO 

10 - 

.1 1 —9 —1 

erg cm s 

1.9 ±0.2 

xl (dof) 



1.29 (395) 


“ Relative abundances from lwilms. Allen & McCi^ j2000l) . 

Gaussian energy fixed at 6.4 keV. 

instead was much redder with no sign of emission lines, suggesting 
an eclipse of the source. We investigated this in more details and 
found that the second spectrum was offset in the spatial direction 
by ~4 pixels compared to the other three. It is likely that the tele¬ 
scope tracking was incorrect and thus another star, 4 away from 
the target, entered in the slit. This star is not recorded in the USNO 
A2 or 2MASS catalogues, but is detected in the WISE images and 
has a similar brightness to Swift J0525.6+2416. This would explain 
the redder spectra and the absence of the emission lines. Thus, we 
conclude that Swift J0525.6+2416 is not an eclipsing system. 

We also searched for radial velocity variations in the emission 
lines. Over the duration of 32 min, Ha was found to vary between 
-163 and -212 km/s (±15 km/s), giving marginal evidence of a 
radial velocity shift. The spectrum shows evidence of interstellar 
reddening with diffuse interstellar bands (DIBs, see Fig. Ill at 6283 
A and possibly at 5780 A. We then used the combined spectrum 
to measure the equivalent width of the DIB feature at 6283 A as 
I.0±0.1 A implying a c olour excess Eb-v in the range 0.39-1.31 
using the relationship of ICordiner et alj ( l201lh . From the hydrogen 
column density of the total absorber derived from the X-ray spectral 
fits we derive E(B-V)= 0.44±0.04 using Nh /E(B-V)= 6.8x10^^ 
( iRvter. Cesarskv & Audouz'3ll975h . We assume the latter one as 
the source colour excess. 

The optical counterpart is also reported in the 2MASS cat¬ 
alogue as 2MASSJ05252270+2413332 with J= 14.65 ± 0.07, 
H= 14.16 ± 0.08, and K= 14.01 ± 0.04 magnitudes. It is also 
detected at longe r wavelengths by th e Wide-field Infrared Survey 
Explorer (WISE: IWright et alj|201(]|) as J052522.84+24133333.6 
with WI (3.4^) = 13.62 ± 0.05 and W2 (4.6^) = 13.49 ± 0.06. 
The 2MASS J-H and H-K colour indexes corrected for intervening 
interstellar absorption are compatible with those of main sequence 
stars w ith spectral types in the range s G8-K3 and G7-M1, respec¬ 
tively iStraizvs & Lazauskaitd 1200^ . The WISE colours (Wl - 
W2)= 0.13±0.08 that can be considered free from intervening ab¬ 
sorption imply a colour temperature of 4000±200 K compatible 
with a donor star in the range K5-M0. However, we must note that 
the accretion disc even in the nIR could mimic a late-type stars 
with different colours. Therefore, only a lower limit to the distance 
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Two new Intermediate Polars 1 
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Figure 4. The spectrum of Swift J0525.6+2416 derived by combining three 
spectra obtained using WHT and ISIS on the night of 2008 Oct 6. The 
spectrum has been calibrated with a standard star to remove the instrumental 
response, but is not on an absolute scale. The Balmer series of Hydrogen ai‘e 
seen in emission as is He II (4686 A). Hq, (6562.8 A) is the stronger line 
in emission. Vertical lines mark the position of diffuse interstellar bands 
(DIBs) at 5780 A(marginal detection) and 6283 A. 
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of about 600 pc can be given, assuming that the 2MASS K band 
absol ute magnitude is that of a Ml star with Mk = 4.93 (iKniggel 

l2006h . 


Swift J0525.6+2416 can be fitted with a multi-temperature opti¬ 
cally thin emission with complex absorption. Also in this source 
we find that the spin pulses are mainly due to photoelectric absorp¬ 
tion originating in the pre-shock accretion flow, which is consistent 
with the accretion-curtain scenario. We found an indication of a 
subsolar abundance of oxygen. We derive a rather massive WD of 
1.0 M 0 . Although most hard X-ray IPs have been found to possess 
WD masses not so different from those of other CVs, this two new 
members appear instead to harbour rather heavy primaries. 
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4 CONCLUSIONS 

IGRJ04571+4527 shows a coherent X-ray signal with a period of 
1222.6 ± 2.7 s, consistent with that found at optical wavelengths, 
confirming it as the WD spin period of the accreting WD. We 
found indication of an orbital modulation of about 8 h, which is 
close to the longer ( 6.2 h) pe riod found in optical spectroscopy 
llThorstensen & HalDernll2()13h . The X-ray spectrum is well fitted 
with a multi-temperature optically thin emission with complex ab¬ 
sorpti on as frequently found in magnetic CVs (e.g. lBernardini et^ 
I 2 OI 2 I) . We conclude that IGRJ04571+4527 is an IP. Spectral 
changes with respect to the spin rotation of the WD indicates that 
the absorption originates from material in the pre-shock flow. The 
spin pulses are consistent with the generally accepted accretion- 
curtain geometry, wh ere material flows towards the W D poles in 
an arc-shaped curtain teosen. Mason & Cordoval|l988h . When the 
primary pole points towards the observer the optical depth of the in¬ 
flowing material above t he shock is larger produ c ing an X-ray flux 
minimum. We note that iThorstensen & Halneml j2013h estimate a 
high interstellar reddening for this source, EfB-Vj-^LO, which im¬ 
plies N_tf ~ 6.8 X 10^^ cm“^. They suggested that this source is 
locate d farther than 500 pc as initially proposed by iMasetti et alj 
( I 2 OI 0 I) . The value we found for Ah is a factor 2 lowe r and s eems 
to favour the initial distance estimate of IMasetti et alj j2010l) . We 
derive a massive W D of 1.12 M 0 , which is at the high end of C V 
WD masses (see e.g. lFerrario, de Martino & Gansickql20151) . 

Swift J0525.6+2416 shows a coherent X-ray signal at 226.28 
s that we identify as the WD spin period. Thus also this source is 
a CV of the IP type. There is no evidence for orbital variability 
in both the X-ray and optical XMM-Newton light curves. From a 
re-analysis of previous optical spectra we do not confirm the pres¬ 
ence of an eclipse. Further optical investigation is encouraged to 
constrain the orbital period. The broad band X-ray spectrum of 
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